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C$ theory taking 2+1 dimensional QED as an example. This theory possesses confinement 



A potential between mesons is extracted from 4-point functions within lattice gauge 



Qh and dynamical fermions. The resulting meson-meson potential has a short-ranged hard 
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repulsive core due to antisymmetrization. The expected dipole-dipole forces lead to at- 
traction at intermediate distances. Sea quarks lead to a softer form of the total potential. 
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1. Introduction. The extraction of an effective interaction, or potential, between two 
composite hadrons from a quantum field theory is of fundamental importance. On one 
hand, since a potential is a central object of classical and quantum mechanics, its definition 
based on a quantum field is interesting in itself. On the other hand, the potential between 
clusters of particles has proved very useful in many studies of subatomic physics. Here 
we have in mind the quark-gluon degrees of freedom at the subnuclear level which are 
thought to give rise to the nucleon-nucleon forces usually described by meson theory. The 
vacuum of QCD contains both virtual gluons and quarks. At short distances the nucleon- 
nucleon forces are mediated by gluon exchange between the constituent quarks whereas 
for longer distances the production of quark-antiquark pairs should be the dominating 
mechanism. Quark-antiquark exchange can be treated as an effective meson exchange 
leading to the construction of the Bonn and Paris potentials [f]. Gluon exchange can be 
studied by phenomenological potential and bag models which provide some insight into 
the interaction mechanism of the six-quark system [2]. 

The aim should be to calculate the nucleon-nucleon forces directly from the quantum 
fields of QCD. Such a project has to deal with the effective forces between composite colour 
singlets, in other words, their "chemistry". This goes along with a characteristic change 
of scale leaving a residual effective interaction I0~ 2 to I0~ 3 times smaller than a typical 
hadron mass. The numerical difficulties involved in going from the GeV to the MeV scale 
in a lattice computation are evident. For this reason it seems appropriate to investigate 
simpler lattice models which exhibit similar physics. Naturally, one needs to insist on 
confinement and the presence of fermions. In this spirit we report on a simulation aimed 
at extracting an effective meson-meson interaction from a QED gauge field model in 2 + I 
dimensions with staggered fermions which has confinement. 

In the last years lattice calculations with static quarks have demonstrated that the 
potential between two three-quark clusters yields an attractive potential limited to the 
overlap region of the baryons [3, 4]. A hard repulsive core of the potential, as stipulated 
by experiments and their interpretation, could not be observed in the region where the 
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two nucleons have relative distance close to zero. Beyond the static theory, employing 
dynamical quark propagators, correct antisymmetrization and quark exchange become 
possible [5]. In this framework the mechanism which leads to a repulsive core will be 
investigated here analytically and numerically within lattice QED 2 +i- 

2. Lattice QED. We have studied a euclidean lattice field theory in 2 + 1 space-time 
dimensions for a U{Y) gauge symmetry group and staggered fermions. The action of the 
system is S(U, x, x) = Sg(U) + Sf(U } x, x) which depends on the gauge fields U and on 
the fermion fields x, X- F° r the gauge action Sg we use Wilson's plaquette formulation 

Sg = P E (1 " ^ R e T r U^x)U v (x + ^(x + v)Ul(x)) , (1) 

where f3 = 2N/g 2 contains the gauge coupling g and the number of colours N which is 
one in our case. For the fermionic action Sf we employ the Kogut-Susskind prescription 
[6] 

S F = |X) r ^ \{x{ x )U ^{x)x{x + p.) - x{x + vWl^xix)) + m F ^2x{x)x{x) \ , 

(2) 

where Nf is the number of external flavours, nip the dynamical quark mass and the Y Xfl 
are reminiscent of the Dirac matrices. 

Expectation values of physical quantities, defined through path integrals, are estimated 
via Monte Carlo simulations. We compare results for QED 2 +i with and without dynamical 
fermions. The simulation was done on an L 2 X T = 24 2 X 32 periodic lattice with the 
compact U(l) Wilson action at f3 = 1.5 and staggered fermions with Nf = and Nf = 2, 
both with nip = 0.1. We used 64 independent gauge field configurations generated with 
the molecular dynamics algorithm [7]. Fermion propagators were computed utilizing a 
random source technique [8] with 32 and 16 random sources in the quenched case and in 
full QED2+1, respectively. 

3. Meson-Meson Correlator. The one-meson field is a product of staggered Grass- 
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mann fields x an d X with external flavours u and d 

<i>x{t) = Xd{xt)xu{xt) . (3) 
The meson-meson fields with relative distance r = y — x are then constructed by 

$At) = L- 2 J2J2 S r,y-^M t ) ( l ) y{ t ) = L- 2 J2xd{xt)xu{xt)xd{r + x t) Xu (r + x t) . (4) 

x y x 

The information about the dynamics of the meson-meson system is contained in the 
4-point time correlation matrix 

cf,(Mo) = <$J(*)M*o)> - <$J(*)X<M*o)> , (5) 

where ( ) denotes the gauge field configuration average. Working out the contractions 
between the Grassmann fields the following diagrammatic representation is obtained 

£(4) = C {4A) + C {4B) _ C {4C) _ C {4D) ^ 

= 1 1 + X - M - M • < 7 » 

Each of the four contributions to the correlator comprises the exchange of gluons and 
sea quarks. For diagrams C^ AA ^ and C^ AB ^ those take place between the mesons, whereas 
diagrams C^ AC ^ and C^ AD ^ correspond also to interaction mediated by the exchange of 
valence quarks. Denoting contractions of the Grassmann fields by 

n n 

...X ...x ■■■ = G n , with n = 1 . . . 4 , (8) 

we have for example 

43 21 12 34 

C (AA) -{^^,^) = {G 1 Cr 2 G,G\) 1 with f=y-x,f' = y'-x', (9) 



where ~ stands for the sums with normalization factors that carry over from (4). The 
gauge configuration average is taken over the product of all four propagators G. 

In order to define an effective meson-meson interaction it is crucial that the nonin- 
teracting components in are isolated. This is achieved by means of a cumulant (or 



cluster) expansion [9] of the gauge field average. For example 

+ (G 1 G* 4 )(G 3 G* 2 ) + ((GrGlGsGD) , (fO) 

where the last term defines the cumulant. In the first term no interactions from gluons or 
sea quarks between the mesons exist. It is thus interpreted to describe freely propagating 
mesons on the lattice. With reference to (8)— (10) define 



21 12 43 34 

C {4A) = <4<<y <4'<M ~ (GiG* 2 )(G 3 Gl) . (11) 



Diagram 

C {4B) 

can be analyzed in a similar fashion, which then leads us to define C^ AB \ 
The sum of those 

(7(4) = + (12) 

respects the boson symmetry on the meson-field level and constitutes the free meson- 
meson correlator with the relative interaction switched off. 

4- Effective Interaction from Time Evolution. The deviation of from C'^ 
contains the residual effective meson-meson interaction. To find its definition, one may 
calculate the correlation matrix C between two elementary boson fields cf>(xt) that live on 
the lattice sites. The corresponding quantum mechanical Hamilton operator consists of a 
free and an interacting part, H = H -\-Hi. It is possible to derive an explicit expression for 
Hi in terms of the correlators of order n = and n = 1 from the perturbative expansion 
of C fulfilling 

C(t - t ) = e -Hi(t-to) . ( 13 ) 

This result leads us, by way of analogy, to define an effective meson-meson interaction as 



Ul = ~Tt 



, with C = C^C^C^ . (14) 

t = t 



This definition is meant as a relation between matrices and holds independently of the 
basis. It can be tested both in momentum space [10] and in coordinate space. The above 



equations build a bridge between the quantum field theoretical correlation functions on 
the lattice and an effective quantum mechanical Hamiltonian. 

5. Hard Core from Adiabatic Approximation. An analysis of the coordinate-space 
matrix elements of the correlators C^ and C^ reveals an interesting perspective on the 
repulsive core. From C^ expressed in terms of the fermion propagators it is easy to see 
analytically that 

c~2(t, t ) = Cff, \t, t ) + Cff, \t, t ) — Cfiti \t, t ) — c\t ? i \t, t ) = o , 

if f = or r ' = . (15) 

In our simulation we choose the relative distance between the mesons to be the same at the 
initial and final times of the propagation, r = r'. In the spirit of the Born-Oppenheimer 
approximation it is assumed that the interaction proceeds faster than the motion between 
the mesons [11]. This seems to be justified because the quark mass rap = 0.1 is relatively 
large. In this sense we replace the sum over all eigenstates in the spectral representation 
of the correlation function by an average term with an effective energy W(f) 

C£\t,t ) = WW)\ 2 e- En{t - to) ~ c(r)e- w ^-^ . (16) 

n 

To the extent that the energy W(f) of the meson-meson system at fixed relative distance 
r can be extracted from the large-t behaviour of the diagonal elements we may conclude 
from (15)— (16) that W(f = 0) = +oo, provided the strength factor c(r ) is nonzero at 
r = 0. In this case the effective interaction possesses a hard repulsive core. It is due to 
the anticommuting nature of the constituent fermion fields leading to Pauli repulsion. 

6. Results. We want to test numerically the two ways suggested above through (13)- 
(14) and (16), respectively, to extract an effective meson-meson potential from a quantum 
field theory. The correlators C'^ and C'^ were computed for time slices t = 15 ... 19 
around the symmetry point t c = 17 in periodic time. The potentials were obtained 
from cosh-fits to the correlators C. The first definition (13)— (14) represents a matrix and 
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Figure 1: Meson-meson potentials as a function of relative distance r for the direct chan- 
nels (A+B) extracted from quenched QED 2 +i- The result (a) 7i-i(r) from the effective 
time evolution is compared with (b) W/(r) of the adiabatic approximation. In both cases 
the direct channels lead to a repulsive potential barier. Curves are guide to the eye and 
error bars result from a jack-knife analysis. 

leads to a nonlocal Hamilton operator 7ii(r } r ') where r, r ' is the separation between 
the mesons. Since computation of the entire 4-point correlation matrix is very time- 
consuming we only obtained the diagonal elements. This also allows a direct comparison 
with W/(r) = W(r) — 2m from the second definition (16), where 2m = 1.031(9) is the 
mass of the noninteracting two-meson system. 

We begin the description of our simulation data with the quenched theory, Np = 0. 
The result for 7i-i(r) including the direct diagrams (A+B) only is presented in fig. fa 
and that for W/(r) is plotted in fig. fb. Both results correspond to a repulsive potential 
barier. The quark-exchange diagrams (C+D) yield similar behaviour (not shown) except 
their contributions enter with a minus sign into the total correlations (6), (7). 

In fig. 2 all four diagrams (A+B)-(C+D) are included. For both 7i-i(r) and W/(r) 
there is some evidence for a repulsive core, as discussed earlier, followed by attraction at 
medium distances. The shapes emerge as a subtle interplay between all graphs (A+B) 
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Figure 2: Full meson-meson potentials as a function of relative distance r for all four 
graphs (A+B)-(C+D) in pure QED 2 +i. Both total potentials (a) 7i-i(r) and (b) W/(r) 
hint at a repulsive core at short range and exhibit attraction at medium distances. The 
resulting potential originates from a sensitive interplay between all contributing graphs. 
Curves are phenomenological fits to the data (see text) with jack-knife error bars. 

and (C+D). We tried a phenomenological fit function of the type 

Vj(r) = or" 7 - f3e~ Sr , (17) 

and found the following parameters: a = 19.7[10~ 3 ], f3 = 55.2[10~ 3 ],7 = 2.35,<5 = 0.254, 
for Hi, and a = 37.2[1CT 3 ], fi = 81.4[1CT 3 ],7 = 5.42,5 = 0.316, for W/. An exploratory 
solution of the radial Schrodinger equation in two dimensions to the corresponding poten- 
tial form indicates a loosely bounded state. In both approximations additivity between 
the potential from the quark-exchange graphs (C+D) and that of the direct term (A+B) 
with respect to the total potential, - ft( c + D ) ~ n (A+B)-(c+D) ^ ig fulfiUed rea _ 

sonably well. This is a nontrivial point since the extraction of the potentials from the 
corresponding correlation functions is a nonlinear procedure. 

Simulations with Np = 2 external flavours of dynamical fermions also have been 
performed. In this case 16 random sources (instead of 32) were used to generate the 
quark propagators. Further, the potentials from the correlator ratios are smaller. For 
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Figure 3: Full meson-meson potentials as a function of relative distance r for all four 
graphs (A+B)-(C+D) extracted from QED 2 +i with two flavours. Both total potentials 
(a) 7i-i(r) and (b) W/(r) resemble those of the pure gauge case. Dynamical fermions 
lead to softer interactions. The precision of the computation is not as accurate as in the 
quenched simulation. 

these reasons the data are considerably more noisy. Both simulations were performed at 
the same value of f3 = 1.5. For ease of comparison, in the following pictures we rescale 
the energy axis to that of the quenched computation by a factor equal to the mass ratio 

2m(N F = 0) 1.031(9) 



1.65. 



18 



2m(N F = 2) 0.63(1) 

In fig. 3a the results for the total potential 7i-i(r) are shown. There seems no qualitative 
change compared to the pure gluonic case, except that sea quarks apparently lead to a 
softer behaviour. The total potential W/(r) from the adiabatic approximation is shown 
in fig. 3b. Again it resembles Tii and dynamical quarks lead to a certain softening of 
the interaction mechanism. Fits to the functional form (17) converged only for Wj (the 
dashed curve has simply been scaled). 



7. Conclusion. We have proposed two possibilities to extract an effective interaction 
between two composite particles. This brings us closer to our aim of extracting poten- 
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tials of two nucleons consisting of three quarks each from lattice QCD. Such a project 
necessitates the computation of a 4-point hadron Green function which is equivalent to 
a correlation function of six quark propagators within the QCD path-integral. A feasible 
approximation, which we study at present, might be to make two quarks of a baryon 
heavy so that they only act as spectator quarks in the scattering process [12]. 

This paper contains a first study of the influence of both gluons and dynamical quarks 
on meson-meson systems in the framework of QED 2 +i- The sea quarks constitute the 
meson exchange in the sense of the Yukawa theory. With the extension of the pure gauge 
action to the full QED action the contribution of virtual mesons was investigated. It 
was found that gluon exchange and valence-quark exchange are equally important. Sea 
quark effects seem to be relatively small but might be decisive for the definite shape of 
the potential. In the model used mesons feel a short-ranged hard-core potential due the 
Pauli exclusion principle acting between their fermionic constituents. The medium-range 
residual forces between colour-neutral quark-clusters are attractive. It should be noted 
that a computation of meson-meson scattering phase shifts based on Liischer's proposal 
[13] indicates the same behaviour of the potential [14]. The correct antisymmetrization 
of the meson-meson correlator was important for the emergence of a hard core. This is in 
contrast to earlier work [3] in QCD with static quarks where antisymmetrization is not 
possible. In the case of SU (3) the nonabelian colour structure will not necessarily lead 
to a vanishing correlator, as for U{Y) in (15), when the relative nucleon-nucleon distance 
approaches zero. Preliminary studies seem to indicate attraction. Further investigations 
with light valence quarks are necessary to solve the hard-core problem of the nucleon- 
nucleon interaction. 
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